Rhombic preordering on a square substrate 
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A competition of incommensurate symmetries occurs whenever a system is forced to conform to an ordering 
that is different from the intrinsically preferred structure of the system itself. As a model system of such a com- 
petition, we study the rivalry between the triangular ordering of hard disks and the square symmetry induced by 
a periodic square substrate. By using density functional theory as well as Monte Carlo computer simulations, we 
determine the full phase behavior for the case of one particle per minimum. We observe a rhombic preordering 
structure preceding the hexagonal solid as a direct consequence of the competing symmetries. Furthermore, the 
square-rhombic transition is reentrant with increasing substrate interaction. Our predictions can be verified in 
experiments of colloids in laser fields. 
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If a monolayer of particles is adsorbed on a substrate with 
an incompatible symmetry, there is a competition between 
the phase that the adsorbed particle would form without sub- 
strate and the phase that reflects the symmetry of the substrate. 
As already described by Landau [1], there must be a first or- 
der transition between phases with incompatible symmetries. 
Phase transitions induced by substrates are widely studied, 
for example on substrates with a one dimensional commen- 
surate structure [2], which enhances the triangular order that 
is present without substrate. On other substrates with more 
complex periodic or even aperiodic symmetry, new phases 
have been observed [3-11]. For small interactions between 
the particles and the substrate, the phase preferred by the par- 
ticles in the absence of a substrate strength prevails, while the 
phase preferred by the substrate is enforced for a strong sub- 
strate strength. In many systems, additional phases at interme- 
diate substrate strength have been observed, for example for 
nanoparticles or micron-sized colloids [3, 6, 9, 10], for vor- 
tices in type-II superconductors [12, 13], for adsorbed atoms 
(cf., e.g., [7, 14, 15]) or molecules (for a recent review, see 
[16]). An intermediate phase also has been predicted from 
defect-mediated melting theory [17]. 

In this letter we study the competition of incompatible sym- 
metries in a minimal model system. We consider hard disks 




FIG. 1 : Schematic view of the square substrate potential of lattice 
constant a. The disks have diameter a and on average there is only 
one particle located at each minimum of the substrate. 



adsorbed on a square substrate for a density that is chosen such 
that there is always one disk per minimum of the substrate. A 
schematic depiction of the model is shown in Fig. 1. Without 
the substrate potential, a system of hard disks at any tempera- 
ture forms a triangular solid (with quasi-long range order) [18] 
if the density is sufficiently high. Evidently, the disks prefer 
triangular order which is incompatible to the square symme- 
try of the substrate. By using a fundamental-measure density 
functional theory (DFT) as well as Monte Carlo simulations, 
we determine the phase behavior depending on the packing 
fraction and the strength of the particle- substrate interactions. 
At low packing fraction (large wave length of the external 
potential), the system is in a modulated square fluid phase. 
When increasing the packing fraction, we observe a rhom- 
bic preordering as a consequence of the incompatibility of the 
symmetries prior to the first order transition into the triangu- 
lar phase at high packing fraction. In the rhombic ordering, 
the particles are still bound to the minima of the square sub- 
strate. However, every second row or every second column 
of particles is shifted in one direction, while the other rows 
or columns are displaced in the opposite direction. In this 
way, the rhombic structure allows for a larger mean distance 
between the particles acting as a precursor to the hexagonal 
solid, where the distance between nearest neighbors is maxi- 
mal. Using the mean- field type DFT, we find a second-order 
square-fluid to rhombic transition, while with Monte Carlo 
simulations there is no first or second order transition but a 
continuous crossover from the modulated square symmetry to 
rhombic ordering. 

Our results indicate a more general scenario in which 
two competing incommensurate structures lead to the emer- 
gence of a preordering structure which inherits properties 
from both competing phases. As a further result, a reentrant 
square-rhombic transition or crossover is found for increasing 
strength of the substrate interaction. Our model system can 
be realized experimentally by sterically-stabilized colloids on 
a flat surface or at an air-water interface exposed to a laser 
interference pattern with square symmetry, similar to previ- 
ous experiments with colloidal particles where laser fields 
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were used to model an external potential [2-4, 6, 10, 19]. 
The limit of strong square substrates has already been stud- 
ied [3, 5, 11, 20], usually with the emphasis on the case of 
multiple colloids per minimum [3, 5, 6, 11]. 

In our model, the hard disks have a diameter cr which serves 
as length scale. The thermal energy ksT is the energy scale. 
The external potential is 



(a) 5 



Kxt(r) = Fo 1 



(1) 



with the amplitude Vb, the four reciprocal lattice vectors 
{kj} = {(±1,±1)^}, and the lattice constant a of the sub- 
strate square lattice. We study the system at unit filling, 
i.e. there is one particle per square, see again Fig. 1. Unit 
filling implies that the areal number density is which 
immediately translates into the dimensionless area fraction 
T] = 7rcr^/4a^ of the system. 

In the following, we shall first present results for the phase 
diagram from fundamental-measure DFT [21-24] in the two- 
dimensional plane spanned by the area packing fraction r] and 
the substrate potential strength Vo/kBT. Subsequently, we 
describe our Monte Carlo simulation results. 

In DFT [25], a grand canonical free energy functional 
(^(T, /i, A, [p(r)]) is minimized with respect to the density 
profile p{r) at fixed temperature T, chemical potential /i and 
given area of the system A. The functional 1^(T, /i, A, 
is conveniently split into three parts, 

n[p{r)] = TiM^)] + ^exc[p(r)] + ^ext[p(r)], (2) 

where Tid[p{r)] = ksT J drp(r) [in {A^p{r)) - l] is the 
exactly known ideal gas functional which includes the (irrele- 
vant) thermal wavelength A. J-'exc[p(r)] denotes the nontrivial 
excess free energy functional resulting from the interaction 
between the particles for which we adopt the recently devel- 
oped fundamental-measure theory for hard disks [24]. The 
third term describes the coupling to the external potential of 
Eq. (1) and reads 



J drp(r)[K. 



t(r)-/i]. 



(3) 



We use a Picard iteration [23] scheme to minimize the grand- 
canonical free energy Q freely on a fine grid. The chemical 
potential p is used as a Lagrangian multiplier so that the av- 
erage number of particles (i.e. the packing fraction r]) in the 
system is fixed. We analyze the symmetry in the minimiz- 
ing equilibrium density profile and obtain the phase diagram 
shown in Fig. 2. For vanishing external potential Vb = 0, 
we recover the bulk hard disk freezing transition from a dis- 
ordered fluid to a hexagonal crystal. The freezing transition 
is first order in this approximation of the DFT [24] . Increas- 
ing the amplitude of the external potential Vb, particles in the 
fluid phase arrange on the square substrate and thus form a 
modulated fluid phase, see the density profile in Fig. 2, where 
black crosses display the potential minima. Consequently, the 
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FIG. 2: (a) Phase diagram for hard disks on a square- substrate show- 
ing the three different phases, square modulated fluid, rhombic, and 
triangular crystal. The dashed line indicates the second order and 
the solid line the first order phase transition, where the dot denotes 
the triple point. Density contour plots of the three phases are in- 
cluded obtained for a fixed external potential with V^lkuT — 1 and 
packing fractions (a') ry = 0.7069, [(a"), (b)] ry = 0.7257 , and (c) 
r\ — 0.7383. The stars indicate the simulation parameters for Fig. 3. 



phase transition shifts to higher packing fractions, as the exter- 
nal potential suppresses the triangular structure. For reduced 
amplitudes higher than Vb/Zc^T = 0.82, a rhombic ordering 
of particles occurs. The rhombic structure can be thought of 
as a shift of even numbered rows of particles away from the 
lattice sites of the square lattice in one of the four possible di- 
rections, while odd numbered rows are shifted in the opposite 
direction. Thereby the square to rhombic transition can be in 
principle continuous, similar to a Martensitic transition [26]. 
In fact, our DFT studies show that the square-to-rhombic fluid 
transition is second order while the crystallization of either 
modulated fluid or rhombic phase into the triangular lattice 
is first order due to the incompatibility of the two structures. 
The resulting density contour plot of the triangular crystal (see 
Fig. 2) shows a distorted hexagonal crystal. The reentrant 
behavior can be understood by considering the interplay be- 
tween entropy and the substrate potential: At low interactions 
with the substrate, the fluid is nearly unperturbed, except for a 
density modulation prescribed by l^xt(r). When the substrate 
interaction is not too strong, the entropy, which is mainly lim- 
ited by interactions with the neighbors, can be increased by 
forming the rhombic structure, which has a larger mean inter- 
patricle distance. At even higher Vb, the fluctuations are sup- 
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pressed to a point where the particles hardly interact with each 
other, such that they simply conform again with the structure 
enforced by the substrate. 

In our Monte Carlo computer simulations, we fixed the 
number of particles N, area A = No? and temperature T. 
Typical configurations for three Vb values, along with the cor- 
responding scattering function or structure factor S'(k) [27] 
are shown in Fig. 3 for a packing fraction r] = 0.69, well 
below the bulk phase transition. For a weak substrate po- 
tential, the system forms small clusters with hexagonal order, 
which are also found in the fluid for Vb =0. The field aligns 
the clusters, such that rows of particles occupy rows of min- 
ima in the substrate potential in the x and y directions. The 
resulting scattering profile has twelve-fold symmetry as pre- 
dicted by Nelson [17], which is the result of a superposition 
of two scattering profiles each exhibiting hexagonal local or- 
der, where the second profile is rotated by 90° compared to 
the first. For Vo/kBT = 2, the structure is completely dif- 
ferent: the competition between the substrate potential, that is 
minimized by a square structure, and the free volume that is 
larger for triangular local structure, results in a locally rhom- 
bic structure. Formation of equilibrium rhombic clusters with 
their rows along the two different directions causes the rhom- 
bic order to be finite ranged. The shifted rows cause large 
peaks at k^^ = (27r/a) (1,1/2) in between the main peaks. 
Finally for large Vo/ksT = 4, the particles are forced to be 
near one of the minima of the substrate potential, which pre- 
cludes rhombic ordering. Comparable results from DFT can 
be found in the Supplemental Material [27] . 

It is clear from the scattering profile that the peak height 
at ki^ is a good order parameter for rhombic ordering. We 
show SCkn) as a function of Vb for different system sizes in 
Fig. 4(a). For every system size considerably larger than the 
largest rhombic cluster size in the system (i.e. N > 32^), the 
scattering function has no finite size dependence, an unam- 
biguous sign of exponential decay of the rhombic positional 
order with distance. In the inset of Fig. 4(a), we show the 
corresponding peak height S\]<.r) as resulting from a Fourier 
analysis of the density field in DFT; for details see [27]. It is 
nonzero in the rhombic phase, because the rhombic positional 
order is long-ranged in DFT. We see that the rhombic order, 
measured by 5'(ki^), in the simulations is much larger around 
the area where the DFT predicts a stable rhombic phase, so 
we would expect to find phase transitions where the S{]<.r) 
changes most sharply in the simulations. In Fig. 4(b), we show 
the average potential energy, which would show a discontinu- 
ity at a first order phase transition, and the mean squared en- 
ergy fluctuations divided by the number of particles TV, which 
would show a peak that increases with the system size at a sec- 
ond order phase transition. Neither quantities show an indica- 
tion of a phase transition, although the variance of the energy 
shows system size-independent peaks around the Vb values, 
where the rhombic order sharply changes. Clearly, we have 
found a very unusual transition, which has a clear structural 
signature, but no thermodynamic footprint. Furthermore, as 
mentioned before, the Monte Carlo simulations confirm the 




FIG. 3: Typical configurations from Monte Carlo simulations for 
a system with 128^ particles (only a small section of the system is 
shown) and the corresponding scattering function (insets) for a pack- 
ing fraction 77 = 0.69 and external field strength (a) Vo/kBT = 1, 
(b) 2 and (c) 4 (as indicated in Fig. 2). The particles are col- 
ored according to their local structure: black denotes no rhom- 
bic order, orange and blue denote perfect rhombic order with rows 
along the x and y directions, respectively. The scattering peak at 
ki? = (27r/a) (1, 1/2) due to the rhombic pre-ordering is indicated 
in the inset of (b). 
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reentrant behavior found in the DFT. 

In conclusion, we have employed DFT calculations to 
study the phase behavior of hard disks on a square substrate 
and found rhombic preordering which was confirmed in our 
Monte Carlo simulations. Therefore, we have shown that DFT 
can be successfully applied to freezing on incommensurate 
structures. Only the order of the phase transition between the 
square to rhombic transition was underestimated by the mean- 
field like DFT. 

Our results show how freezing occurs on a substrate with a 
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structure that is incompatible with the symmetry of the bulk 
crystal. We found a new intermediate ordering that is a com- 
promise between the two symmetries. The model that we 
studied applies to the competition between square and hexag- 
onal structures where the new emerging intermediate phase 
is rhombic, but we anticipate that this is a general scenario. 
Hence, we expect a comparable phase behavior for other com- 
peting symmetries, as shown for a rectangular substrate in 
the Supplemental Material [27]. However, it is not immedi- 
ately self-evident how a possible preordering, which corre- 
sponds to the rhombic phase on the square lattice, might look 
like on substrates with other symmetries. Therefore, an in- 
teresting question arising from our results is whether the in- 
termediate phases reported for substrates with incompatible 
symmetries, e.g., the phases with 20-fold bond orientational 
order [9] or the Archimedean-like tiling phases [10] on sub- 
strates with decagonal quasicrystalline symmetry, are remi- 
niscent of a preordering before a first order transition between 
completely incompatible symmetries. Furthermore, it will be 
interesting to explore other model systems where the incom- 
patibility of different symmetries is induced by the bound- 
aries [28-30] or is due to an incommensurate filling fraction 
of the substrate minima (see, e.g. [3, 6]), and to study the 
modification of friction due to the different phases [31]. 

We thank M. Oettel for helpful discussions. This work was 
supported by the DFG via SPP 1296 and SFB TR6. 
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FIG. 4: Monte Carlo simulation results for packing fraction 77 = 0.69 
and for system sizes N = 32^, 64^, 96^ and 128^. (a) The scattering 
function 5'(k) at the k vector that is indicative of rhombic order, 
= 27r/a ( 1 , 1 / 2) , as a function of the amplitude Vb of the external 
field. S'^kn) obtained from DFT in the inset, (b) The fluctuations 
of the energy divided by the system size as a function of Vb- The 
inset shows the average energy per particle, for which the error is of 
order of the line thickness. 
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Simulation details 

The scattering function S'(k) (or structure factor) shown in the insets of Fig. 3 of the main text is calculated 
using a fast Fourier transform Pinst(k) of the instantaneous density profile of the corresponding snapshot 
binned on a two-dimensional grid: 5'inst(k) = |pinst(k)p/A^. The method of plotting the scattering profile 
results in a slight smoothing of the peaks. As a result, the peaks at k = (27r/a) (±1, ±1), which are actually 
essentially point-like, obtain a finite width rendering them visible. For the results in Fig. 4(a) of the main 
text, a more precise calculation of the average scattering function was used: 



5(kfl) 



N 



(1) 



which was averaged over multiple configurations during the simulation and over the different symmetry 
related k^^ vectors. The coloring of the particles in Fig. 3 of the main text is performed as follows, the degree 
to which particle i are orange and light blue is equal to Sx^\i) and Sy'^\i) repectively, where 



max -— 

a=±l TV. 



COS 



(2) 



and S'y'^'ii} is given by the same expression with x and y interchanged. The sum over j is limited to the Ni 
particles that are in a square of width 5a centered around the position of particle i. 

Direct comparison of the results from simulations and theory 

In order to clarify the similarities of simulations and theory, first the density profiles obtained from DFT are 
computed for the state points marked by black stars in Fig. 2(a) and the results are shown in Fig. SI. Similar 
to the simulation snapshots shown in Fig. 3, particles in the square phase at high interaction strength Vq are 
strongly pinned or in other words, the density peaks are very high (see Fig. Sl(c)). At low Vq (Fig. SI (a)), 
the square phase is a modulated fluid with much lower density peaks. In DFT, the rhombic phase spans the 
whole system (see Fig. Sl(b)), whereas a simulation snapshot consists of regions of rhombic and regions of 
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FIG. SI: Density contour plots of 6x6 squares for fixed packing fraction 77 = 0.69 and external field strength (a) 
Vo/kBT — 1, (b) Vo/kBT — 2 , and (c) Vb/ZcsT = 4, similar to those indicated by stars in Fig. 2. 
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FIG. S2: Structure factor 5'(k) obtained from DFT at the state points indicated in Fig. SI 

square symmetry. As the state point for the rhombic phase is very close to the phase boundary, the rhombic 
structure is hardly visible but it can be confirmed by analyzing the profiles. 

In analogy to the analysis of the results from our MC simulations, a structure factor related quantity can 
be calculated. The Fourier transform of the mean density distribution p(k) is used as 

^'(k) = (3) 

conatining a normalization factor l/5'o. Here, the average and the Fourier transform are taken in opposite 
order compared with Eq. 1. The results are in good agreement with those presented in Fig. 3. There is an 
intermediate peak at k^^ = ^(1/2, 1) only in the rhombic phase. In the simulations, the shift of particles 
from the square to rhombic phase occurs in the x as well as the y direction forming different clusters. In 
contrast, DFT yields only one way of dislocation. Thus, the intermediate peaks can only occur in one of the 
two possible directions. In the inset of Fig. 4(a), we plot the peak height at k^^ as a function of Vq/ZcbT. All 
S''(k) are normalized with a factor l/5'o such that the peak height obtained from DFT agrees with that of 
the black curve in Fig. 4(a) resulting from our simulations. 

Additional to the analysis of a square substrate, an external potential with rectangular symmetry can be 
induced for a system with on average one particle per substrate minimum. This setup also yields a rhombic 
preordering as shown in a density contour plot in Fig S3. 




FIG. S3: Rhombic phase at packing fraction 77 = 0.69 on a rectangular pattern with size ratio 1:1.2 and external field 
strength Vo/kBT — 1.5 indicated by black crosses. 



